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SECTION  1 


THE  MULTIPLE  PHASE-SCREEN  METHOD 


In  the  Multiple  Phase-Screen  method  (MPS),  the  ionized  medium  is  divided 
into  a  finite  number  of  layers  and  each  layer  is  represented  by  a  thin  diffracting  screen 
(see  Appendix  A  and  References  3  and  4  for  more  details).  MPS  provides  an  exact 
solution  to  the  parabo'ic  wave  equation  describing  a  plane,  unmodulated  carrier-wave 
traversing  this  series  of  phase-perturbing  screens.  Hence  MPS  is  “exact”  in  that  the 
fidelity  of  its  results  is  limited  only  by  the  extent  to  which  the  parabolic  equation 
approximation  and  phase-spectrum  description  of  the  environment  preserve  the  true 
dynamical  and  environmental  characteristics  of  the  problem. 

It  should  be  stressed  that  MPS  calculates  actual  values  of  signal  amplitude 
and  phase  (I  and  Q  voltage  levels)  at  the  receiver  plane;  the  method  is  fundamentally 
different  from  procedures  that  replace  the  ionized  layer  with  an  explicit  statistical 
model  for  the  received  signal.  In  MPS,  although  realizations  of  the  phase-changing 
screens  are  generated  in  a  random  manner  based  upon  a  sequence  of  pseudo-random 
numbers  (consistent  with  the  layer’s  phase-variance  description),  once  these  random 
numbers  are  generated,  the  phase-screens  are  completely  defined  and  known.  Simi¬ 
larly,  the  wave  field  propagated  through  the  phase-screens  is  subject  to  exact  calcula¬ 
tion.  However,  a  different  sequence  of  pseudo-random  numbers  will  yield  a  different 
set  of  phase-screen  realizations  and  hence  a  different  realization  of  the  received  electric 
field  Realizations  of  the  received  signal  amplitude  and  phase  are  important  results  of 
the  simulation  and  in  turn  serve  as  direct  input  to  detailed  receiver  simulations. 

MPS  can  handle  all  levels  of  ionospheric  disturbances,  from  the  least  severe 
with  only  minor  phase  fluctuations  to  the  most  severe  cases  of  frequency  selective 
Rayleigh  fading.  It  is  in  the  cases  of  weak  scattering  that  MPS  methods  prove  most 
useful,  since  in  these  cases  the  signal  amplitude  statistics  are  neither  Rayleigh  nor 
Rician  -  their  description  is  quite  complex  and  cannot  be  accurately  generated  by 
ot  her,  more  analytically-based,  models  that  make  assumptions  regarding  the  statistical 
character  of  the  signals. 

Though  exact,  the  MPS  method  is  limited  in  its  application  by  its  require¬ 
ment  that  the  electric  field  and  phase  be  adequately  represented  by  values  at  a  discrete 
number  of  grid  points;  that  is,  in  order  to  provide  a  faithful  description  of  the  phase 
distribution  and  avoid  angular  aliasing  and/or  scattering,  the  phase-screen  grid  point 
spacing  must  be  fine  enough  and  the  phase-screen  length  must  be  long  enough.  For 


1 


propagation  environments  parameterized  by  large  values  of  phases variance  and/or 
small  fluctuation  scale-size,  these  requirements  may  call  for  a  prohibitively  large  num¬ 
ber  of  grid  points. 

Fortunately,  there  exist  other  methods  for  simulating  strongly-scattered  sig¬ 
nals.  The  Antenna/Channel  Impulse  Response  Function  (ACIRF)  code  (see  Refer¬ 
ence  l),  which  implements  the  DNA-approved  statistical  generation  technique,  is  far 
less  computationally  intensive  than  MPS  but  relies  strongly  on  the  assumption  that 
Rayleigh  statistics  govern  the  amplitude  fluctuations  of  the  signals  of  interest. 

Thus  the  MPS  method  is  well-suited  for  the  cases  of  comparatively  weak 
scattering,  where  the  Rayleigh  assumptions  upon  which  ACIRF  is  based  break  down, 
and  ACIRF  provides  appropriate  realizations  for  strongly-scattered  Rayleigh  signals, 
when  the  grid  resolution  requirements  of  MPS  prove  too  costly  computationally. 
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SECTION  2 


THE  DATA  BASE 


2.1  WHAT  IS  REPRESENTED  IN  THE  DATA  BASE. 


The  magnetic  tape  contains  the  results  of  a  series  of  fifteen  MPS  runs. 
Each  run  or  “case”  corresponds  to  a  specific  choice  for  electron  density  variance  and 
radiofrequency  wavelength,  and  the  output  of  each  MPS  run  is  a  series  of  signal 
realizations  -  actual  values  of  signal  amplitude  and  phase  (I  and  Q  voltage  levels)  -  at 
the  receiver  plane.  The  simulated  cases  fall  conceptually  into  three  distinct  wavelength 
classes,  each  with  its  own  unique  set  of  five  fluctuation  levels.  Within  each  wavelength 
class,  the  electron  density  variance  values  specified  produce  scintillation  effects  that 
range  from  weak  to  strong  scattering: 


Electron  Density  Variance:  cm  6] 

Wavelength 

3  cm 

9.5  cm 

30  cm 

2.12  x  1013 

4.26  x  10u 

9.37  x  109 

6.72  x  1013 

1.35  x  1012 

2.96  x  1010 

2.12  x  1014 

4.26  x  1012 

9.37  x  1010 

6.72  x  1014 

1.35  x  1013 

2.96  x  10u 

2.12  x  1015 

4.26  x  1013 

9.37  x  1011 

The  ionized  layer  through  which  the  plane  wave  propagates  imparts  a  phase- 
shift  that  is  a  function  of  distance  along  the  horizontal  axis  of  the  layer.  The  variance 
of  this  phase-shift,  the  “phase-variance” ,  is  determined  in  part  by  the  electron  density 
variance,  fluctuation  scale  sizes,  and  the  layer-thickness.  The  MPS  method  replaces 
the  entire  layer  by  a  smaller  number  of  thin  layers,  and  compresses  the  thin  layers 
into  phase-perturbing  screens  that  p.eserve  phase- variance. 

The  index  of  refraction  fluctuations  of  the  ionized  layer  are  characterized  by 
a  two-component  power  spectral  density  (see  Appendix  A)  whose  parameters  include: 
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Outer  Scale  Size  -  10  km/rad, 

Freezing  Scale  Size  -  300  m/rad, 

Inner  Scale  Size  -  2  m/rad, 

Intermediate  Scale  Parameter  n  -  1.75,  and 
Transition  Scale  Parameter  it'  -  2.5. 


Together  with 


Layer  Thickness  for  an  Individual  Phase  Screen  -  24.6  km, 
Total  Number  of  Layers  -  11,  and 
Briggs  Parkins  Angle  -  90  deg., 


the  power  spectral  density  of  the  layer  (and  electron  density  variance)  determines  the 
power  spectral  density  of  the  MPS  phase-screens. 

The  MPS  simulation  is  also  dependent  upon  geometry.  Relevant  descriptors 
are  presented  in  Figure  1.  There  are  217  grid  points  along  the  space  axis  for  a  spatial 
resolution  of  two  meters. 

MPS  requires  as  input  the  phase-variance  per  screen  and  the  wave  frequency. 
The  preceding  table  becomes  (see  Appendix  A  for  conversion  details): 


Phase- Variance  per  Screen:  [rad2] 

Frequency 

10 

GHz 

3.16 

GHz 

1  GHz 

5.63 

X 

101 

1.13 

X 

101 

2.49 

1.78 

X 

102 

3,59 

X 

101 

7.85 

5.63 

X 

102 

1.13 

X 

102 

2.49  x  101 

1.78 

X 

103 

3.59 

X 

102 

7.85  x  101 

5.63 

X 

103 

1.13 

X 

10s 

2.49  x  102 

Table  1  summarizes  the  resulting  MPS  calculations.  The  entries  for  each 
case  are  average  values  over  ten  realizations.  For  each  frequency  class,  as  the  total 
phase-variance  increases,  the  mean-squared  value  of  the  logarithm  of  signal  amplitude, 
(x2),  increases  from  .010  toward  .5,  the  normalized  intensity  variance,  S4,  becomes 
partially-saturated,  and  the  decorrelation  distance,  £o>  decreases. 

Figures  2  and  3  are  typical  realization  profiles.  In  Figure  2,  a  weakly- 
scattered  case,  intensity  fluctuations  are  due  to  a  mixture  of  refraction  and  diffraction 
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Ax 
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Figure  1.  MPS  geometry. 
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Table  1.  Summary  of  MPS  calculations. 


°N « 

Electron 

Density 

Case  # 

Frequency 

Variance 

[cm-6] 

4001 

10  GHz 

2.12E+13 

4002 

(A  =  3  cm) 

6.72E+13 

4003 

2.12E+14 

4004 

6.72E+14 

4005 

2.12E+15 

4006 

3.16  GHz 

4.26E+11 

4007 

(A  =  9.5  cm) 

1.35E+12 

4008 

4.26E+12 

4009 

1.35E+13 

4010 

4.26E  +  13 

4011 

1  GHz 

9.37E+09 

4012 

(A  =  30  cm) 

2.96E+10 

4013 

9.37E+10 

4014 

2.96E+11 

4015 

9.37E-I-11 

Total 

Phase 


Variance 

[rad2] 

<  x2  > 

5. 

fo 

[m] 

6.20E+02 

.00980 

.201 

191 

1.96E+03 

.0312 

.369 

105 

6.20E+03 

.103 

.698 

58.0 

1.96E+04 

.315 

1.12 

32.4 

6.20E-F04 

.570 

1.31 

18.1 

1.25E+02 

.00983 

.201 

457 

3.95E+02 

.0315 

.366 

243 

1.25E+03 

.107 

.678 

132 

3.95E+03 

.327 

1.08 

73.1 

1.25E+04 

.559 

1.24 

40.7 

2.73E+01 

.00986 

.200 

1120 

8.64E+01 

.0316 

.361 

563 

2.73E+02 

.109 

.648 

296 

8.64E+02 

.331 

1.01 

160 

2.7E+03 

.548 

1.18 

88.3 
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CASE  4002  (x2  *=  .03)  realization  1 


Distance  (Km) 


20- 

m 

33  io. 


IVIRC 


o. 


C/1 

£  -io. 


-20, 


10. 


20.  30.  40.  50. 

Distance  (Km) 


60. 


70. 
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Figure  2.  Plots  of  phase,  intensity,  and  their  respective  power  spectral 
densities  for  a  realization  with  (x2)  =  .03  (10  GHz). 


CASE  4004  (x  ~  -32)  realization  1 


Distance  (Km) 


and  consist  of  small  spikes  with  magnitude  less  than  5  dB.  Signal  phase  as  a  function 
of  distance  is  smooth  compared  to  Figure  3,  which  is  a  more  strongly  scattered  case. 

In  Figure  3,  the  phase  fluctuations  are  much  more  extreme  and  there  is 
a  slow  drift  up  from  0  to  30  km  and  then  down  again.  This  “TEC-wandering”  or 
“phase-wandering”  is  due  to  large  scale  structure  in  the  intervening  layer. 

The  intensity  plot  of  Figure  3  displays  frequent  10  dB  spikes  and,  more 
importantly,  some  instances  of  extremely  sharp  fades.  For  the  strongest-scattering 
case  simulated,  these  intensity  fluctuations,  due  primarily  to  diffraction,  begin  to 
approach  Rayleigh  as  can  be  seen  by  comparing  the  plots  of  Figure  4. 

Figure  5  is  a  graph  of  average  phase  power  spectral  density  at  the  receiver 
corresponding  to  the  most  weakly-scattered  case  at  10  GHz.  It  is  apparent  that 
the  spectral  density’s  slope  is  governed  by  2n'-l  for  large  spatial  wavenumber  and  by 
2n—  1  for  smaller  K,  where  n  and  n'  are  the  respective  intermediate  and  transition  scale 
parameters  that  help  define  the  index  of  refraction  fluctuation  power  spectral  density. 
The  small  “hump”  at  K  =  3  x  10~2  is  due  to  a  purely  geometrical  Fresnel-filtering 
effect;  it  corresponds  to  the  first  Fresnel  zone. 

Evidence  that  the  grid-spacing  was  chosen  small  enough  to  preserve  phase 
information  is  given  by  the  relatively  small  number  of  times  that  the  phase  changed 
by  more  than  90  degrees  from  grid-point  to  grid-point.  For  the  worst  case,  the  case 
of  highest  frequency  and  largest  phase  variance,  there  were  only  242  of  these  “phase- 
slips”  out  of  more  than  138,000  points. 


2.2  HOW  TO  ACCESS  INFORMATION  IN  THE  DATA  BASE. 


The  MPS  calculations  produced  150  output  files;  that  is,  ten  realizations  for 
each  of  the  fifteen  cases  of  frequency  and  phase-variance.  Analysis  of  these  files  re¬ 
vealed  little  variation  between  realizations  within  a  given  case.  Therefore,  considering 
the  size  of  these  files,  and  in  keeping  with  previous  procedures,  only  a  subset  of  the 
total  150-file  set  was  placed  in  the  data  base.  The  realizations  chosen  for  inclusion 
were  based  on  the  calculated  (x2)  value  for  the  case,  where  (x2)  is  the  mean-square 
value  of  the  logarithm  of  signal  amplitude.  Included  for  each  case  are  the  following 
files: 

1.  MPS(case_number).AVG  -  the  realization  in  which  the  measured  (x2)  fell  closest 
to  the  average, 
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Avg  Phase  PSD 


10.  20.  30.  40.  50.  60.  70. 

Distance  (Km) 

Theoretical  Rayleigh  signal 


Distance  (Km) 

Figure  4.  Strongly-scattered  case  compared  to  a  theoretical  Rayleigh  < 

signal. 


Figure  5.  Theoretical  slope  structure  in  the  phase  PSD  corresponding  to 
the  most  weakly-scattered  case  (10  GHz). 


2.  MPS(case_number).M AX  -  the  realization  with  largest  measured  (x2),  and 

3.  MPS(case_number).Ml\  the  realization  with  smallest  measured  (x2). 

Tables  2a-2c  summarize  the  data  generated.  Each  filename  appears  along 
with  its  layer’s  total  phase  variance  and  measured  values  of  (x2),  £4  index,  and  decor¬ 
relation  distance. 

The  output  files  themselves  each  contain  a  simulation  identification  record, 
records  containing  problem  definition  data,  and  floating  point  realization  data  records, 
which  are  written  as  buffered  blocks  of  (I,Q)-vo!tage  level  pairs.  Appendix  B  provides 
details  concerning  the  output  file  format  and  Appendix  C  lists  a  Fortran  subroutine 
that  may  be  used  to  read  realization  data  into  a  simulation  code. 


Table  2a.  Description  and  summary  of  10  GHz  data  base. 


Carrier  Frequency  =  10  GHz 


Filename 

(x2) 

to 

[rad2] 

H 

average 

6.20E+02 

.00980 

.201 

191 

mps4001.avg 

.00987 

.201 

193 

mps4001.max 

.0104 

.206 

197 

mps4001.min 

.00940 

.195 

190 

average 

1.96E+03 

.0312 

.369 

105 

mps4002.avg 

.0314 

.370 

106 

mps4002.max 

.0330 

.378 

107 

mps4002.min 

.0302 

.355 

104 

average 

6.20E+03 

.103 

.698 

58.0 

mps4003.avg 

.104 

.695 

53.4 

mps4003.max 

.110 

.722 

59.4 

mps4003.min 

.0991 

.675 

57.9 

average 

1.96E+04 

.315 

1.12 

32.4 

mps4004.avg 

.315 

1.12 

32.3 

mps4004.max 

.327 

1.16 

33.1 

mps4004.min 

.298 

1.09 

31.3 

average 

6.20E+04 

.570 

1.31 

18.1 

mps4005.avg 

.569 

1.29 

18.9 

mps4005.max 

.578 

1.34 

16.7 

mps4005.min 

.558 

1.27 

17.6 

12 


Table  2b.  Description  and  summary  of  3.16  GHz  data  base. 


Carrier 

Frequency 

=  3.16  GHz 

Filename 

(x2> 

[rad2J 

M 

average 

1.25E+02 

.00983 

.201 

457 

mps4006.avg 

.00983 

.200 

495 

mps4006.max 

.0104 

.206 

411 

mps4006.min 

.00933 

.195 

453 

average 

3.95E+02 

.0315 

.366 

243 

mps4007.avg 

.0316 

.363 

259 

rnps4007.max 

.0333 

.375 

250 

mps4007.min 

.0299 

.356 

242 

average 

1.25E+03 

.107 

.678 

132 

mps4008.avg 

.107 

.664 

140 

mps4008.max 

.115 

.708 

136 

mps4008.min 

.0989 

.659 

132 

average 

3.95E+03 

.327 

1.08 

73.1 

mps4009.avg 

.327 

1.03 

77.3 

mps4009.max 

.346 

1.11 

73.6 

mps4009.min 

.313 

1.07 

71.3 

average 

1.25E+04 

.559 

1.24 

40.7 

mps4010.avg 

.561 

1.30 

40.7 

mps4010.max 

.597 

1.29 

37.6 

mps4010.min 

.533 

1.20 

42.5 
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Table  2c.  Description  and  summary  of  1  GHz  data  base. 


Carrier  Frequency  =  1  GHz 


Filename 

<x2> 

s4 

to 

[rad2] 

M 

average 

2.73E+01 

.00986 

.200 

1120 

mps4011.avg 

.00987 

.199 

1180 

mps4011.max 

.0108 

.210 

989 

mps4011.min 

.00902 

.190 

1210 

average 

8.64E+01 

.0316 

.361 

563 

mps4012.avg 

.0317 

.357 

573 

mps4012.max 

.0351 

.379 

504 

mps4012.min 

.0290 

.339 

607 

average 

2.73E+02 

.109 

.648 

296 

mps4013.avg 

.109 

.655 

301 

mps4013.max 

.121 

.671 

269 

mps4013.min 

.101 

.598 

317 

average 

8.64E+02 

.331 

1.01 

160 

mps4014.avg 

.324 

.990 

156 

mps4014.max 

.372 

1.06 

147 

mps4014.min 

.306 

.969 

155 

average 

2.73E+03 

.548 

1.18 

88.3 

mps4015.avg 

.545 

1.18 

85.4 

mps4015.max 

.605 

1.27 

81.2 

mps4015.min 

.515 

1.13 

86.4 
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APPENDIX  A 


THEORETICAL  OVERVIEW 


The  scalar  wave  equation, 


(V2  +  k2n2)E  =  0 


where 


V2 

k 


d 2  ,  d2 
dx 2  dz2 
2n 

T  ’ 


and  n  —  1  -f  An  is  the  index  of  refraction  with  An  a  small  perturbation,  describes 
a  plane  wave  propagating  in  the  z-direction  through  a  medium  assumed  to  consist  of 
random  striations  that  are  infinitely  extended  in  the  y-direction  (out  of  the  page  in 
Figure  6a). 

This  assumption,  which  effectively  reduces  the  analysis  to  a  two-dimensional 
geometry,  is  often  a  good  one  since  the  striated  ionization  tends  to  align  itself  with  the 
earth’s  magnetic  field  and  propagation  before  the  Rayleigh  limit  is  reached  is  primarily 
two-dimensional  (the  simulation  can  be  easily  extended  to  handle  a  three-dimensional 
geometry  but  would  require  considerably  greater  computational  effort). 

In  Figure  6b,  the  ionized  region  has  been  partitioned  into  a  number  of  thin 
layers,  each  perpendicular  to  the  direction  of  propagation.  Consider  the  ith  layer, 
centered  at  z  =  z,,  with  thickness  Az. 

In  order  to  analyze,  first,  the  effect  of  the  layer  upon  the  wave’s  phase,  the 
electric  field  is  written  in  terms  of  complex  amplitude  and  phase: 

E(x,z,u>)  =  U(x,  z,ui)e~lkz 


This  has  the  effect  of  separating  the  free  space  propagation  portion  of  the 
electric  field  from  the  phase-changing  portion  that  is  affected  by  the  medium.  Substi¬ 
tuting  into  the  scalar  wave  equation  and  making  order  of  magnitude  arguments  leads 
to  the  parabolic  wave  equation: 

tJ  dLJ 

— 5  -  2ik— —  +  2k2  An(i,z,w)t/  =  0 
dx  dz 
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Suppose  the  irregularities  in  each  thin  layer  are  now  compressed  into  in¬ 
finitesimally  thin  screens  that  preserve  the  phase-perturbing  effect.  The  ionized 
medium  is  thus  represented  by  a  series  of  thin  diffracting  phase-screens  separated 
by  free  space  (Figure  6c). 


For  small  Az,  the  equation  for  propagation  through  the  thin  layer  is  obtained 
from  the  parabolic  wave  equation  by  neglecting  the  first  term;  that  is,  propagation 
effects  within  the  thin  layer  have  been  ignored.  The  analysis  focuses  only  on  those 
terms  governing  the  phase  change  induced  by  the  medium’s  varying  index  of  refraction. 
The  resulting  equation  is  easily  solved  by  separation  of  variables: 


U+{x,Zi,uj) 


LL(x,  Zi,  u>)  exp 


An(x,z', 


(1) 


where  £/_  and  U+  denote  the  values  of  U  as  it  enters  and  exits  the  phase-screen  at 
Zi.  Equation  1  states  that  the  exiting  wave’s  phase  is  a  simple  sum  of  phase  change 
associated  with  the  layer  and  the  phase  of  the  electric  field  at  the  entrance  to  the 
layer. 


Turning  now  to  the  free-space  propagation  portion  of  the  analysis,  the  last 
term  of  the  parabolic  wave  equation  can  be  ignored,  and  the  resulting  equation  is 
easily  solved  by  Fourier  transform  techniques: 

U_(K,zi+l,uj)  =  U+(Kizi,u)ex p  [iir2(zi+1  -  z,)/2fc] 

where  the  circumflex  denotes  the  Fourier  transform  of  U . 

Taking  the  inverse  transform  leads  to 

U-(xtZi+uu)  =  J  U+{K,  z,-,w)  exp  \iK2(zi+l  -  z,)/2&  +  x/fx]  d.K 

or  in  terms  of  the  electric  field 

E-(x,Zi+u u)  =  J  E+{K,Zi,u>)exp[iK2(zi+1  -  z,)/2A:  +  i/fi]  d,K  (2) 

Propagation  of  the  electromagnetic  wave  from  Z{  to  zi+1  has  been  decom¬ 
posed  into  two  parts:  first,  a  pure  phase  change  contribution,  induced  by  irregularities 
in  zC s  layer,  calculated  by  adding  the  phase  change  indicated  by  Equation  1  to  the 
phase  of  the  electric  field  at  the  entrance  to  the  phase-screen  associated  with  z,  and 
second,  a  free-space  propagation  contribution  given  by  Equation  2,  which  is  carried 
out  by  Fourier  transforming  the  phase-perturbed  field,  multiplying  in  K- space  by 
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e*K3(zi+i -zi)/2k,  an(j  pourjer  transforming.  This  two-step  procedure  is  repeated  at 
each  screen,  with  the  final  free-space  propagation  terminating  at  the  receiver  plane. 

Specific  realizations  of  the  random  phase  in  Equation  1, 


4>{x)  =  —k  An(x,z')dz' 


that  define  a  phase-screen  are  obtained  by  sampling  a  distribution  of  phase  shifts 
whose  statistical  properties  are  specified  by  the  power  spectral  density  of  the  electron 
density  irregularities.  Specifically,  the  one-dimensional  phase  PSD  is  related  to  the 
three-dimensional  in-situ  electron  density  fluctuation  PSD  as  follows: 

First,  the  three-dimensional  electron  density  fluctuation  PSD  is  given  by 
(see  Reference  5  for  complete  definitions  of  terms) 

pen  (K  K  ...  8^!(A^)AT3(n,„',ii)iri,i1rW  , 

fbUNc{Kx,Kv,Kz)  -  - -  ffn  ~  3/2) -  / 

[(1  +  KlL\  +  K*L*  +  K\L\  +  2LxyKxKy  +  2LXZKXKZ  +  2LyzKyKz)n 
(1  +  K\t\  +  KX  +  K\t\  +  2  txvKxKy  +  2ixzKxKz  +  2  tyzKyKz)n'-n] 

where  normalization  terms  have  been  introduced  to  make 


I 


dK 


PSDn—z=AN? 


By  making  use  of  the  Markov  approximation  (i.e.,  by  assuming  that  PSDNt 
is  delta-correlated  in  the  z-direction),  the  two-dimensional  differential  phase  spectrum 
can  be  directly  related  to  the  three-dimensional  electron  density  fluctuation  PSD: 

^(Kz,Ky)  =  r7X7PSDNt(Kx,Ky,Kz  =  0) 


Finally,  the  one-dimensional  phase  PSD  can  be  obtained  by  carrying  out 
the  integration  over  the  spatial  wavenumber  Ky  (again,  Reference  5  contains  complete 
definitions  of  terms): 


dP*(Kx)  =  f°°  dP*(K  Ky)dKy 

az  J_00  dz  2tt 

do l _ 2y/frN2T{n  -  l/2)(l  +  6.4/n')Lz. 

T(n  -  1)(1  +  6.4/n)(l  +  L2K7z)n-5{d7  +  t7K7z)n'-n 
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A  practical  technique  for  generating  the  (discrete)  Fourier  transform  of  a 
given  phase-screen  realization  makes  use  of  the  discrete  spectral  density  of  the  desired 
phase: 


4>(mAk)  =  rm\S(mAk)L/27r\1^2 


where  rm  =  ^/l/2(<7Im  +  ig2rn ),  dim  and  g2m  are  samples  from  a  standard  normal 
distribution,  and  S(mAk)  is  the  discrete  spectral  density  corresponding  to  a  given 
layer;  i.e., 

dP+ 


A  Pa  = 


dz 


A  z 


and  S(mAk)  is  the  discrete  version  of  A 

The  “phase-variance-per-screen”  quantity  input  to  MPS  is  obtained  by  in- 

j  2 

tegrating  differential  phase  variance,-^,  over  the  thickness  of  the  layer  represented 
by  the  screen. 


MPS  can  handle  the  regions  in  propagation  space  that  include  weak  scat¬ 
tering  to  near  Rayleigh  scattering.  Its  use  of  FFTs  make  MPS  relatively  fast  and 
inexpensive,  so  many  cases  can  be  run.  Although  MPS  is  practical  only  for  two- 
dimensional  problems  due  to  computer  considerations,  this  is  not  a  serious  problem 
for  reasons  mentioned  earlier. 

The  most  serious  limitation  on  the  MPS  method  derives  from  the  practical 
limitation  on  the  smallness  of  Ax,  the  transverse  mesh-spacing.  As  the  propagation 
environment  degrades,  the  distance  over  which  significant  phase  or  amplitude  fluctua¬ 
tions  occur  becomes  smaller.  In  practice  it  is  required  that  Ax  <  £0  .  This  restriction 
limits  the  severity  of  effects  that  can  be  accurately  simulated  by  the  MPS  method. 

Fortunately,  calculations  reveal  that  UHF  signals  become  Rayleigh  before 
Ax  becomes  so  small  as  to  make  MPS  calculations  impractical.  In  the  Rayleigh  regime, 
other  simulation  methods,  such  as  ACIRF  (Reference  1),  can  be  employed. 
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APPENDIX  B 


MPS  OUTPUT  FILE  DESCRIPTION 


This  appendix  describes  the  format  of  signal  realization  data  base  files.  For 
a  specific  request,  an  ANSI  standard  tape  will  be  written  with  data  in  ASCI  format 
(or  some  other  mutually  agreeable  format).  This  standard  data  format  was  chosen 
over  a  machine-specific  format  to  alleviate  the  need  for  support  of  the  wide  range  of 
computer  systems  in  use  by  government  contractors.  In  addition  to  the  tape,  a  small 
Fortran  program  will  be  included  to  convert  the  specific  format  to  the  unformatted 
style  used  in  this  appendix.  The  overall  file  structure  is  as  described  below  (note  that 
each  output  file  contains  a  single  realization): 


Record 

1  character  identification  record 

2  floating  point  record  A  (problem  definition  data) 

3  floating  point  record  B  (detailed  problem  definition  data) 

4  record  A  (for  check  purposes) 

5  floating  point  realization  data  record  Cl 

6  record  A 

7  record  C2 
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On  the  following  pages  a  detailed  description  of  each  of  the  above  records  is 
provided.  For  each  record,  the  Fortran  write  statements  used  to  produce  the  record 
are  shown  and  followed  by  a  tabular  description  of  the  variables  written. 
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character  identification  record 


C 

C 


Variable 


NID 

IDENT 


CHARACTER  IDENT *80 
INTEGER  NID 


NID  =  80 

WRITE  (1)  NID , IDENT(i :NID) 


Description 


number  of  characters  in  IDENT  string  (nominally  80) 
alphanumeric  identification  of  simulation 
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floating  point  record  A  (problem  definition  data) 


C 

C 

C 


INTEGER  I.N1 

REAL  RDAT1 

DIMENSION  RDAT1 (28) 


N1  =  28 

WRITE  (1)  N1 . (RDATl(I) ,1  =  1  ,N1) 


< 

Variable  Description 


N 1 

RDAT1(1:N1) 

where 

RDATl(  1)  =  SOURCE 

RDAT1(  2)  =  CASEID 
RDATl(  3)  =  FC 

*  RDATI(  4)  =  0. 

*  RDAT1(  5)  =  0. 
RDATl(  6)  =  DLX 

*  RDAT1(  7)  =  0. 
RDATl(  8)  =  CHIBAR 


number  of  floating  point  numbers  in  RDATl  array  (cur¬ 
rently  28) 

floating  point  problem  definition  data 

simulation  source  code  identifier: 

1.  =  MPS 
case  identifier 
carrier  frequency  (Hz) 


decorrelation  distance  in  x-direction  (m) 

(x2)  -  mean  square  log  amplitude  fluctuation  of 
realization 
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Variable 


Description 


RDAT1(  9)  =  S4 
RD  AT  1(10)  =  XLENG 

RDATl(ll)  =  NGRIDX 
RDAT1(12)  =  DELTAX 

*  RDAT1(13)  =  0. 

*  RDAT1(14)  =  0. 

*  RDAT1(15)  =  0. 
RDAT1(16)  =  RNRLX 
RDATl(17)  =  NFREQ 
RDATl(18)  =  FREQO 
RD  AT  1(19)  =  DELTAF 

*  RD  AT  1(20)  =  0. 

*  RDAT1(21)  =  0. 

*  RD  AT  1(22)  =  0. 
RDAT1(23)  =  RANI 
RDAT1(24)  =  RAN2 
RDAT1(25)  =  MAXBUF 


*  RDATl(26)  =  0. 

*  RDATl(27)  =  0. 
RDAT1(28)  =  VERS 


scintillation  index  of  realization 
grid  length  (m) 

number  of  grid  points 
grid  point  spacing  (m) 


number  of  grid  points  per  DLX 
number  of  frequencies 
lowest  frequency  (Hz) 
delta  frequency  (Hz) 


random  number  seed  1 
random  number  seed  2 

maximum  buffer  size  (real  words)  [buffer  size  required 
to  read  largest  realization  record  in  this  file] 


code  version  number 


*  denotes  quantities  not  calculated  by  the  MPS  code.  Realization  data  generated  by 
both  the  MPS  and  ACIRF  computer  codes  share  a  common  file  format.  When  a 
particular  parameter  is  not  calculated  by  the  source  code  the  value  is  assigned  zero. 
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floating  point  record  B  (detailed  problem  definition  data) 


INTEGER  I,N2 

REAL  RDAT2 

C 

DIMENSION  RDAT2(250) 

C 

C 

N2  =  260 

WRITE  (1)  N2 , (RDAT2(I) , 1=1 ,N2) 

I 


Variable  Description 


N2  number  of  floating  point  numbers  in  RDAT2  array  (not 

to  exceed^SO) 

RDAT2(1:N2)  floating  point  detailed  problem  definition  data  (Note 
that  this  data  is  primarily  meant  to  aid  in  detailed  prob¬ 
lem  description  for  internal  MRC  use  and  archive.) 

where  RDAT2(l:N2)  are  yet  to  be  determined  by  the  MPS  and  ACIRF  code  custodi 
ans. 
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floating  point  realization  data  records  (buffered  due  to  size) 


INTEGER  I . K . LFTOVR . MAXBUF ,  N1 , NGRIDX 

REAL  RDAT1 , REALI , REALQ 

C 

PARAMETER  (NGRIDX  =  "number  of  grid  points") 

PARAMETER  (MAXBUF  =  "data  buffer  size  --  must  be  even") 

C 

DIMENSION  RDAT1 (28) .REALI (NGRIDX) .REALq (NGRIDX) 

C 

C 

N1  =  28 

NUMBLK  =  NGRIDX*2/MAXBUF 
C 

DO  I  =  1,  NUMBLK 

WRITE  (1)  N1 , (RDATl(K) ,K=1 ,N1) 

WRITE  (1)  MAXBUF. 

it  (REALI ((I-1)*MAXBUF+K) , REALQ ( (1-1) +MAXBUF+K) ,K=1 .MAXBUF/2) 

ENDDO 

C 

LFTOVR  =  MOD  (MAXBUF ,NGRIDX*2) 

IF  (LFTOVR. GT.O)  THEN 
WRITE  (1)  N1,(RDAT1(K) ,K=1,N1) 

WRITE  (1)  LFTOVR, 

it  (REALI (NUMBLK*MAXBUF+K) , REALQ ( NUMB LK*MAXBUF+K) ,K=1 .LFTOVR/2) 

ENDIF 
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Variable 


Description 


N 1  number  of  floating  point  numbers  in  RD  ATI  array  (cur¬ 

rently  28) 

RDATl(l:Nl)  floating  point  problem  definition  data  (each  of  these 
records  is  the  same  as  the  2nd  record  in  the  file  and 
is  used  as  a  means  of  verifying  that  the  data  file  has 
been  read  in  correctly) 

NGRIDX  number  of  grid  points 

MAXBUF  buffer  length  o"  realization  data  record  (words) 

NUMBLK  number  of  realization  data  record  blocks  (those  of  length 

MAXBUF) 

LFTOVR  number  of  left  over  realization  data  (only  written  when 

MAXBUF  not  a  factor  of  NGRIDX) 

Note  that  the  following  data  were  written  as  I,Q  pairs 
(II, Ql,  12, Q2,  13, Q3,  ...) 

REALI()  array  of  floating  point  realization  data  (I  voltage) 

REALQ()  array  of  floating  point  realization  data  (Q  voltage) 
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APPENDIX  C 


FORTRAN  SUBROUTINE  FOR  SAMPLING  REALIZATIONS 


This  appendix  lists  a  sample  Fortran  subroutine  for  use  in  reading  MPS  or 
ACIRF  realization  data  into  a  simulation  code.  The  prospective  user  may  find  the 
following  advisory  comments  useful. 

If  Ax  is  the  reciever  plane  mesh  spacing,  then  there  are  mesh  points 
per  decorrelation  distance  (note  that  Lq  and  Ax  are  located  in  positions  6  and  12 
respectively  of  the  array  RDATl;  see  the  description  of  floating  point  record  A  in 
the  previous  appendix).  For  a  given  decorrelation  time  r0,  the  mesh  point  index 

corresponding  to  a  particular  time  t  of  interest  is  given  by  ITIME  = 

The  subroutine  interpolates  I,Q  data  from  values  at  bracketing  mesh  points. 

If  the  sample  desired  falls  beyond  the  extent  of  the  realization  data,  then 
the  subroutine  wraps  the  realization  around  to  simulate  a  longer  run. 

If  one  prefers  to  access  the  realization  data  in  the  spatial  domain  rather 
than  in  the  time  domain,  one  can  simply  call  SAMPLE  with  the  time  arguement 

set  appropriately;  i.».,  t  =  (&Ax).  Here,  the  quantity  (fcAx)  is  distance  in  the 

receiver  plane,  and  k  is  the  mesh  point  index. 


SUBROUTINE  SAMPLE  (IUNIT , FILEIN .RESET, TAUO, TIME, VI ,VQ) 

* 

************************************************************************ 


*  * 

*  Fortran  subroutine  to  get  I  and  Q  voltages  from  MPS  or  ACIRF  * 

*  data  files.  * 

*  * 

*  INPUTS  from  argument  list:  * 

*  IUNIT  =  logical  unit  number  for  input  file  * 

*  FILEIN  =  input  file  name  * 

*  RESET  =  logical  reset  flag  (should  be  set  to  TRUE  by  calling  * 

*  routine  when  a  new  file  is  specified  or  a  new  run  ip  * 

*  to  be  started:  will  be  set  to  FALSE  by  this  routine.)  * 

*  TAUO  =  scintillation  decorrelation  time  (s)  * 
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*  TIME  =  time  of  I-Q  data  relative  to  beginning  of  file  (s)  * 

*  * 

*  OUTPUTS  to  argument  list:  * 

*  VI, VQ  =  inphase  and  quadrature  voltage  samples  * 

*  * 


CHARACTER  FILEIN* (*) , IDENT*80 

* 

LOGICAL  RESET 

* 

DATA  NIQ/O/ 

*  maximum  buffer  size 
PARAMETER  (IBUF  =  4096) 

* 

DIMENSION  RDAT1 (28) . RDAT2 (250) . REALI (0 : IBUF/2) . REALQ (0 : IBUF/2) , 
b  CKDAT(28) 

* 

SAVE  RDAT2 . REALI . REALQ . DELTAT , IMIN . IMAX . N1 . NGRIDX . NTOTAL 

* 

************************************************************************ 

* 

*  check  reset  flag  (new  file  or  new  case) 

1  IF  (RESET)  THEN 

CLOSE  (UNIT=IUNIT) 

OPEN  (UNIT=IUNIT ,FILE=FILEIN ,STATUS= ’ OLD  * ,FORM= ’UNFORMATTED ’ ) 

*  read  ID  record 

READ  (IUNIT)  NID,IDENT(1 :MIN0(NID,80)) 

*  read  record  A 

READ  (IUNIT)  N1 , (RDAT1 (I) ,1=1 ,MIN0(N1 ,28)) 

*  read  record  B 

READ  (IUNIT)  N2 , (RDAT2 (I) ,1=1 ,MIN0(N2 ,250) ) 

IF  (NID.GT.80)  WRITE  (*,*)  ’ID  record  exceeds  80  characters’ 

IF  (N1.GT.28)  WRITE  (*,*)  ’Record  A  exceeds  28  words’ 

IF  (N2.GT.250)  WRITE  (*,*)  ’Record  B  exceeds  250  words’ 

IF  (NID.GT.80  .OR.  N1.GT.28  .OR.  N2.GT.250)  STOP 
MAXBUF  =  INT(RDAT1 (25) ) 

IF  (MAXBUF. GT. IBUF)  THEN 

WRITE  (*.*)  ’Buffer  size  must  be  increased  in  SAMPLE’ 

WRITE  (*.*)  ’Set  IBUF  =  ’. MAXBUF. ’  and  recompile’ 
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STOP 
END  IF 

ELZERO  =  RDAT1(  6) 

NGRIDX  =  RDATl(ll) 

NTOTAL  =  NGRIDX 
DELTAX  =  RDAT1(12) 

DELTAT  =  TAUO*DELTAX/ELZERO 
IMIN  =  0 
IMAX  =  0 

IF  (NIQ.EQ.O)  NIQ  =  IBUF/2 
RESET  =  .  FALSE . 

END  IF 

ITIME  =  INTCTIME/DELTAT)  ♦  1 
IF  (ITIME.LT. 1)  THEN 

WRITE  (*.*)  ’Negative  time  cannot  be  used  in  SAMPLE’ 
STOP 
END  IF 

IF  (ITIME.LT. IMIN)  THEN 

WRITE  (*,*)  'Warning:  Time  out  of  sequence  in  SAMPLE’ 
WRITE  (*,*)  ’Will  assume  that  you  want  to  start  over’ 
RESET  =  .TRUE. 

GO  TO  1 
END  IF 


*  Data  is  in  buffer 

2  IF  (ITIME.LT. IMAX)  GO  TO  4 
IF  (ITIME. GE. NTOTAL)  GO  TO  3 


read  check  record 

READ  (IUNIT,END=99)  NCK . (CKDAT(I) . 1=1 .NCR) 

IF  (NCK.NE.N1)  THEN 

WRITE  (*,*)  ’Incorrect  number  of  check  words’ 
STOP 
END  IF 

DO  I  =  1.  NCK 

IF  (CKDAT(I) . NE . RDAT1 (I) )  THEN 
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WRITE  (*,*)  ’Check  word  error’ 
STOP 
END  IF 
END  DO 


*  Read  I.Q  data  into  buffer 
REALI(O)  =  REALI(Niq) 

REALQ (0)  =  REALQ(NIQ) 

READ  (IUNIT ,END=99)  NN . (REALI (I) , REALQ (I) . 1=1 . NN/2) 

NIQ  =  NN/2 

IMIN  =  IMAX 

IMAX  =  IMAX  NIQ 

GO  TO  2 

* 

*  Wrap  realization  around  for  longer  run 
3  REWIND  (IUNIT) 

NTOTAL  =  NTOTAL  +  NGRIDX 
READ  (IUNIT) 

READ  (IUNIT) 

READ  (IUNIT) 

GO  TO  2 


************************************:>.*********•».**  lit********************;** 


*  Interpolate  I.Q  data  to  specified  time 
4  DT  =  AMOD (TIME , DELTAT)  /  DELTAT 
I  =  ITIME  -  IMIN 

VI  =  REALI(I) * (1 . -DT)  ♦  REALI(I+1)*(DT) 

VQ  =  REALQ (I) * ( 1 . -DT)  +  REALQ(I+ 1) * (DT) 

RETURN 

99  STOP  ’**  ERROR:  END-OF-FILE  ON  READ  IN  SUBROUTINE  SAMPLE  **’ 
END 
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APPENDIX  D 


DISCUSSION  AND  DETAILED  SUMMARY  OF  THE  MPS  RUNS 


ways. 


This  appendix  summarizes  the  results  of  the  MPS  calculations  in  several 


First,  Tables  3a,  3b,  and  3c  (corresponding  to  realizations  obtained  at 
10,  3.16,  and  1  GHz,  respectively)  give  the  measured  (x2)  and  Lq  for  all  realizations 
and  their  theoretical  values.  The  “theoretical”  values  referred  to  were  calculated 
using  methods  implemented  in  PRPSIM  (see  Reference  2).  The  set  of  analysis  points 
along  a  given  line  of  sight  were  defined  by  the  line  of  sight’s  intersection  with  the 
MPS  thin  layers.  Theoretical  values  for  (x2)  were  derived  using  the  weak  scintillation 
calculations  implemented  in  PRPSIM,  while  the  theoretical  Lq  was  found  by  searching 
for  the  |  point  of  the  autocorrelation  function  at  the  receiver. 


r 


The  tables  are  followed  by  graphs  (figure  7)  of  the  results  versus  a %  for 
the  three  frequencies.  The  average  value  for  a  set  of  realizations  obtained  for  a  single 
frequency  and  phase  variance  is  indicated  by  the  centered  “X”  and  the  individual 
realizations  are  indicated  by  the  small  dots.  It  can  be  seen  that  the  agreement  between 
the  theoretical  and  measured  values  of  lo  is  excellent.  For  (x2)  the  agreement  is 
excellent  for  values  less  than  about  .3.  Since  the  theoretical  formula  has  been  derived 
assuming  weak  scintillation,  the  deviation  for  large  (x2)  is  understandable. 


Next,  figures  8-12  represent  realizations  obtained  for  a  frequency  of  10 
GHz  and  are  ordered,  as  are  the  other  sets  of  figures  that  follow,  from  weakest  to 
strongest  scintillation  (smallest  to  largest  phase  variance).  The  figures  indicate  the 
level  of  amplitude  and  phase  variation  across  the  receiver;  they  contain  plots  of  signal 
phase  and  intensity,  and  their  respective  power  spectral  densities,  for  the  realizations 
whose  measured  (x2)  fell  closest  to  the  case  average.  These  realizations  correspond 
to  the  realizations  on  the  tape  labeled  by  MPS(case_number) .AVG.  In  a  similar  fash¬ 
ion,  figures  13-17  represent  realizations  obtained  at  3.16  GHz,  and  figures  18-22 
represent  the  1  GHz  case. 
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For  realizations  with  small  (x2)  (.01  and  .03)  the  fading  and  flaring  intensity 
is  generally  confined  to  within  5  dB  of  the  mean  intensity.  For  stronger  scintillation 
((x2)  >  -10)  the  intensity  is  more  severe;1  as  the  scintillation  begins  to  saturate 
((x2)  ~  .5)  worst  case  fading  and  flaring  is  obtained.  As  the  scintillation  increases 
not  only  does  the  range  of  fading  and  flaring  increase  but  the  separation  of  the  sharp 
fades  or  flares  decreases.  For  the  weak  scintillation  cases  the  phase  PSD  plot  shows  a 
slight  kink  near  10~2  rad/m.  This  kink  represents  the  effect  of  Fresnel  filtering  at  the 
first  Fresnel  zone. 


Finally,  figures  23-27  provide  a  statistical  comparison  of  the  received  power 
obtained  from  the  MPS  realizations  with  three  theoretical  distributions:  (l)  Rayleigh; 
(2)  Rician;  and  (3)  Nakagami-m.  These  comparisons  make  use  of  the  MPS  realizations 
from  cases  4001-4005  (figures  8-12),  whose  carrier  frequency  is  10  GHz.  The  first  plot 
in  each  figure  compares  the  probability  density  functions  (PDF)  for  the  power  fades 
and  flares,  and  the  second  plot  compares  the  cumulative  distribution  functions  (CDF). 

A  histogram  was  accumulated  for  each  of  the  ten  realizations  of  a  given  case. 
In  the  figures,  dots  indicate  the  average  bin  values  over  the  ten  realizations,  while  bars 
show  the  standard  deviation  of  bin- values.  It  can  be  seen  from  the  figures  that  both 
Rician  and  Nakagami-m  distributions  provide  somewhat  better  fits  for  S4  <  1.0  than 
the  Rayleigh  distribution,  with  the  Nakagami-m  distribution  offering  a  slightly  more 
accurate  model  for  the  probability  of  fades. 

For  cases  with  >  1.0,  the  Rician  model  provides  a  better  estimate  for 
the  probability  of  fades  than  the  Nakagami-m,  while  the  Nakagami-m  closely  matches 
the  distribution  of  flares.  The  PDF  shown  in  the  last  figure,  CASE  4005,  reflects 
the  difference  between  the  two  intensity  plots  shown  in  figure  4  (page  10),  where  the 
theoretical  Rayleigh  signal  lacks  the  flares  of  intensity  greater  than  10  dB  that  the 
strongly-scattered  case  possesses.  The  PDF  shows  a  corresponding  difference  between 
the  distribution  of  power  greater  than  5  dB  for  the  MPS  realizations  (dots)  and  the 
Rayleigh  distribution  (solid  line). 

Current  DNA  convention  regards  signals  possessing  (x2)  values  greater  than  0.1  as  being  “Rayleigh” 
signals.  Rayleigh  signals  are  generated  by  propagation  paths  penetrating  significant  nuclear  effects 
regions.  Typically  a  slight  perturbation  (towards  increasing  effects)  in  link/effects  region  geometry  for 
weak,  near-Rayleigh  signals  gives  rise  to  rapidly  increasing  values  for  (x2)-  The  rapid  change  in  (x2) 
implies  a  corresponding  rapid  transition  to  Rayleigh  statistics.  Thus  the  0.1  value  is  a  conservative  yet 
practical  threshold. 
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The  Rayleigh  PDF  and  CDF  are  given  by  the  formulas: 


f(P)  =  exp  (-P)  ,  F(P)  =  1  -  exp  (-P) 


The  Rayleigh  functions  are  indicated  on  the  plots  by  the  short  dashed  lines. 
The  Rician  PDF  and  CDF  are  given  by  the  formulas: 


HP)  - 

F(P)  =  [P  f(P')dP' 

Jo 

(  \[Pq  y/p\ 

=  l-(?  V>  a\ 


Q(a,b)  =  exp{-2!-^}£Ji0(g)k/k(<.6) 


The  appropriate  values  P0  and  o  are  obtained  by  deriving  them  from  the  measured 
scintillation  index 

(p1)  -  (py 
*  ■  (py 

The  formulas  used  are  the  following: 


The  Rician  functions  are  indicated  on  the  plots  by  the  solid  lines. 

The  Nakagami-m  PDF  and  CDF  are  given  by  the  formulas: 
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where  7  is  the  incomplete  Gamma  function  and  m  is  chosen  by  comparing  the  theo¬ 
retical  and  measured  values  of  the  scintillation  index: 


m  = 


si 


The  Nakagmai-m  functions  are  indicated  on  the  plots  by  the  long  dashed  lines. 
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Table  3a.  Detailed  MPS  realization  summary  of  cases  4001-4005  (10 
GHz). 


Total 

Phase 

Case  Variance 
[rad2] 

4001  6.20E+02 


4002  1.96E+03 


4003  6.20E+03 


Realization  (x2)  S4 


[mj 

1 

.00990 

.201 

194 

2 

.00965 

.200 

196 

3 

.00993 

.203 

174 

4 

.00987  (AVG) 

.201 

193 

5 

.00970 

.203 

183 

6 

.00960 

.198 

203 

7 

.00940  (MIN) 

.195 

190 

8 

.00989 

.202 

185 

9 

.0104  (MAX) 

.206 

197 

10 

.00966 

.199 

199 

Average 

.00980 

.201 

191 

Theory 

.00959 

196 

1 

.0315 

.367 

106 

2 

.0307 

.368 

107 

3 

.0315 

.376 

96.0 

4 

.0314  (AVG) 

.370 

106 

5 

.0308 

.376 

101 

6 

.0305 

.364 

111 

7 

.0302  (MIN) 

.355 

104 

8 

.0316 

.372 

102 

9 

.0330  (MAX) 

.378 

107 

10 

.0308 

.365 

109 

Average 

.0312 

.369 

105 

Theory 

.0303 

108 

1 

.104 

.689 

58.6 

2 

.101 

.700 

59.2 

3 

.104  (AVG) 

.695 

53.4 

4 

.105 

.707 

58.5 

5 

.102 

.707 

56.1 

6 

.100 

.687 

61.2 

7 

.0991  (MIN) 

.674 

57.9 

8 

.105 

.704 

56.6 

9 

.110  (MAX) 

.722 

59.4 

10 

.101 

.689 

60.4 

Average 

.103 

.698 

58.0 

Theory 

.0959 

59.7 
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Table  3a.  Detailed  MPS  realization  summary  of  cases  4001-4005  (10 
GHz)  (continued). 


Total 

Phase 

Variance 

Realization 

(x2) 

s4 

^0 

[rad2] 

M 

1.96E+04 

1 

.319 

1.15 

32.7 

2 

.320 

1.11 

33.0 

3 

.317 

1.14 

29.8 

4 

.322 

1.11 

32.6 

5 

.298  (MIN) 

1.09 

31.3 

6 

.306 

1.12 

34.1 

7 

.315  (AVG) 

1.12 

32.3 

8 

.313 

1.10 

31.6 

9 

.327  (MAX) 

1.16 

33.1 

10 

.309 

1.11 

33.7 

Average 

.315 

1.12 

32.4 

Theory 

.303 

33.4 

6.20E+04 

1 

.574 

1.37 

18.3 

2 

.563 

1.28 

18.4 

3 

.578  (MAX) 

1.34 

16.7 

4 

.572 

1.32 

18.3 

5 

.558  (MIN) 

1.27 

17.6 

6 

.573 

1.27 

19.1 

7 

.571 

1.33 

18.1 

8 

.570 

1.31 

17.8 

9 

.568 

1.29 

18.5 

10 

.569  (AVG) 

1.29 

18.9 

Average 

.570 

1.31 

18.1 

Theory 

.959 

18.7 

Table  3b.  Detailed  MPS  realization  summary  of  cases  4006-4010  (3.16 

GHz). 


Total 

Phase 

Case  Variance 
(rad2) 

4006  1.25E+02 


4007  3.95E+02 


4008  1.25E+03 


Realization  (x2) 


(mj 

1 

.00990 

.200 

469 

2 

.00969 

.200 

468 

3 

.0104  (MAX) 

.206 

411 

4 

.00978 

.200 

463 

5 

.00972 

.202 

436 

6 

.00983  (AVG) 

.200 

495 

7 

.00933  (MIN) 

.195 

453 

8 

.00987 

.201 

440 

9 

.0103 

.205 

473 

10 

.00950 

.197 

478 

Average 

.00983 

.201 

457 

Theory 

.0960 

471 

1 

.0318 

.364 

246 

2 

.0311 

.364 

24« 

3 

.0332 

.376 

221 

t 

.0314 

.366 

245 

•J 

.0309 

.372 

232 

0 

.0316  (AVG) 

.363 

259 

7 

.0299  (MIN) 

.356 

242 

8 

.0317 

.368 

235 

') 

.0333  (MAX) 

.375 

250 

10 

.0304 

.357 

253 

Average 

0315 

.366 

243 

Theory 

.0304 

249 

1 

.109 

.671 

134 

> 

.103 

.674 

135 

3 

112 

.695 

121 

1 

108 

.680 

133 

5 

103 

.680 

128 

6 

.107  (AVG) 

.664 

140 

7 

0989  (MIN) 

.659 

132 

A 

108 

.692 

129 

9 

!  1  5  !  M  A  X) 

.708 

136 

lM 

lot 

.659 

138 

A  vei  age 

107 

.678 

132 

I’lieoi  v 

0900 

136 

•11 


Table  3b.  Detailed  MPS  realization  summary  of  cases  4006-4010  (3.16 
GHz)  (continued). 


Total 

Phase 

Variance 

Realization 

<x2> 

to 

[rad2| 

[m] 

3.95E+03 

1 

.346  (MAX) 

1.11 

73.6 

2 

.315 

1.05 

74.5 

3 

.344 

1.10 

67.1 

4 

.331 

1.10 

73.5 

5 

.313 

1.04 

70.6 

6 

.327  (AVG) 

1.03 

77.3 

7 

.314 

1.08 

72.9 

8 

.313  (MIN) 

1.07 

71.3 

9 

.343 

1.09 

74.7 

10 

.320 

1.08 

76.1 

Average 

.327 

1.08 

73.1 

Theory 

.304 

75.1 

1.25E+04 

1 

.571 

1.28 

41.0 

2 

.540 

1.21 

41.4 

3 

.597  (MAX) 

1.29 

37.6 

4 

.569 

1.24 

40.9 

5 

.551 

1.23 

39.4 

6 

.557 

1.23 

43.0 

7 

.561  (AVG) 

1.30 

40.7 

8 

.553 

1.24 

39.9 

9 

.563 

1.20 

41.5 

10 

.533  (MIN) 

1.20 

42.5 

Average 

.559 

1.24 

40.7 

Theory 

.960 

41.9 

42 


Table  3c.  Detailed  MPS  realization  summary  of  cases  4011-4015  (1  GHz) 


Case 


4011 


4012 


4013 


Total 

Phase 

Variance 

Realization 

<x2) 

54 

^0 

[rad2] 

M 

2.73E+01 

1 

.00987(AVG) 

.199 

1180 

2 

.00990 

.198 

1140 

3 

.0108  (MAX) 

.210 

989 

4 

.00981 

.201 

1170 

5 

.0101 

.205 

1050 

6 

.00902  (MIN) 

.190 

1210 

7 

.00974 

.199 

1130 

8 

.00948 

.197 

1060 

9 

.0102 

.204 

1200 

10 

.00969 

.198 

1170 

Average 

.00986 

.200 

1120 

Theory 

.00957 

1190 

8.64E+01 

1 

.0317 

.358 

581 

2 

.0317  (AVG) 

.357 

573 

3 

.0351  (MAX) 

.379 

504 

4 

.0315 

.363 

567 

5 

.0323 

.369 

534 

6 

.0290  (MIN) 

.339 

607 

7 

.0313 

.357 

560 

8 

.0303 

.359 

540 

9 

.0325 

.369 

594 

10 

.0312 

.356 

589 

Average 

.0316 

.361 

563 

Theory 

.0302 

582 

2.73E+02 

1 

.109  (AVG) 

.655 

301 

2 

.107 

.637 

304 

3 

.121  (MAX) 

.671 

269 

4 

.110 

.657 

298 

5 

.110 

.649 

282 

6 

.101  (MIN) 

.598 

317 

7 

.107 

.640 

295 

8 

.104 

.660 

286 

9 

.111 

.667 

307 

10 

.106 

.642 

309 

Average 

.109 

.648 

296 

Theory 

.0957 

304 

43 


Table  3c.  Detailed  MPS  realization  summary  of  cases  4011-4015  (1  GHz) 
(continued). 


Total 

Phase 

Case  Variance 
[rad2) 


4014  8.64E+02 


4015  2.73E+03 


Realization 

<*2> 

to 

lm] 

1 

.351 

1.05 

162 

2 

.317 

.996 

164 

3 

.372  (MAX) 

1.06 

147 

4 

.323 

1.02 

161 

5 

.306  (MIN) 

.969 

155 

6 

.314 

.944 

170 

7 

.340 

1.01 

160 

8 

.324  (AVG) 

.990 

156 

9 

.321 

1.02 

165 

10 

.341 

1.02 

168 

Average 

.331 

1.01 

160 

Theory 

.302 

165 

1 

.554 

1.20 

88.7 

2 

.516 

1.14 

89.5 

3 

.605  (MAX) 

1.27 

81.2 

4 

.538 

1.15 

88.5 

5 

.545  (AVG) 

1.18 

85.4 

6 

.580 

1.17 

93.7 

7 

.565 

1.27 

88.2 

8 

.515  (MIN) 

1.13 

86.4 

9 

.534 

1.12 

90.0 

10 

.522 

1.17 

92.4 

Average 

.548 

1.18 

88.3 

Theory 

.957 

90.7 

44 
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Figure  9.  Representative  phase,  intensity,  and  power  spectral  density 
plots  for  case  4002. 


Intensity  PSD 


Prob  [P/P0  ^  X]  PDF  of  Power 


Figure  23.  Comparison  of  Rayleigh,  Rician,  and  Nakagami-m  PDFs  and 
CDF8  for  case  4001. 
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Prob  [P/P0  ^  X]  PDF  of  Power 


MPS  CASE  4004 


Figure  26.  Comparison  of  Rayleigh,  Rician,  and  Nakagami-m  PDFs  and 
CDFs  for  case  4004. 
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